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We recently studied the formation and evaporation of primordial black holes in a simple brane- 
world cosmology, namely Randall-Sundrum Type II. Here we study the effect of accretion from the 
cosmological background onto the black holes after formation. While it is generally believed that 
in the standard cosmology such accretion is of negligible importance, we find that during the high- 
energy regime of braneworld cosmology accretion can be the dominant effect and lead to a mass 
increase of potentially orders of magnitude. However, unfortunately the growth is exponentially sen- 
sitive to the accretion efficiency, which cannot be determined accurately. Since accretion becomes 
unimportant once the high-energy regime is over, it does not affect any constraints expressed at the 
time of black hole evaporation, but it can change the interpretation of those constraints in terms of 
early Universe formation rates. 



PACS numbers: 98.80.Cq 

I. INTRODUCTION 

Primordial black holes (PBHs) are relics which may 
form in the early Universe [[l] and survive to have effects 
in later epochs, thus shedding light on the physical pro- 
cesses operating during the Universe's early stages. In a 
recent paper (Guedens, Clancy and Liddle ||), we con- 
sidered the physics of black holes in a popular variant to 
the standard cosmology, the Randall-Sundrum Type II 
braneworld cosmology || . We showed that the presence 
of the fifth (AdS) dimension could significantly modify 
the properties of black holes, and indeed that black holes 
evaporating at key epochs of the Universe might behave 
as higher-dimensional objects. 

The results of that paper assumed that after forma- 
tion the dominant process affecting the mass of the black 
holes was Hawking evaporation. In the standard cosmol- 
ogy it is believed that the neglect of accretion of ma- 
terial from the cosmological background is a very good 
approximation, except perhaps immediately after forma- 
tion, though the calculations are rather uncertain. How- 
ever, we remarked that it was not necessarily true that 
accretion would also be negligible during the high-energy 
regime. The purpose of this paper is to investigate accre- 
tion during the high-energy regime. 1 We will find that 
accretion can indeed be the dominant effect in the high- 
energy regime, though accurate calculation of its effects 
is presently impossible. 

Our calculations follow the notation of Ref. 0, and 
we will consider both the standard cosmology and the 
Randall-Sundrum Type II cosmology, which features a 



As we were completing this paper, a paper appeared by Majum- 
dar [W , using the formalism of our earlier paper fcj to investigate 
this same question. Our results are in broad agreement, though 
we will highlight the key differences. 
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high-energy regime where the Friedmann equation is 
modified. The AdS radius of the extra dimension will 
be denoted I, and M4 and I4 indicate the usual four- 
dimensional Planck mass and length (we set c = I). 
The four-dimensional cosmological constant set to zero. 
We use M for the black hole mass and r$ for the 
Schwarzschild radius. If ro < I then the black hole will 
be effectively five-dimensional and the appropriate five- 
dimensional Schwarzschild solution is used. 



II. ACCRETION FORMALISM 

Calculations of PBH accretion in the standard cosmol- 
ogy have a long history but are plagued with signifi- 
cant uncertainties. First of all, while it is known that 
black holes forming during radiation domination cannot 
be much smaller than the horizon size (otherwise pressure 
forces would prevent collapse), it is unclear how much 
smaller they might be. Secondly, it is unclear precisely 
how efficient a black hole might be at absorbing incom- 
ing radiation from a cosmological background, or how 
important the effect of backreaction might be. Results 
prove extremely sensitive to assumptions made concern- 
ing these quantities. While early work by Zel'dovich and 
Novikov J?! speculated that PBHs might even be able to 
grow as fast as the horizon, subsequent work, especially 
by Carr and Hawking Q , made a convincing case that 
such growth could not occur, and moreover that once 
the PBH became significantly smaller than the horizon 
accretion would become very inefficient. This view is now 
widely held, though there are papers where accretion is 
found to be important, including interesting work by Ha- 
cyan Q in which incoming radiation was described by 
a Vaidya metric matched to a flat radiation-dominated 
Universe, and in which the mass was found to grow pro- 
portional to the horizon size (M(t) — ► 0.06 — 0.08 1) at 
late times even for very small initial mass under the ide- 
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alization of perfect radial inflow. 

There are different approaches to estimating the accre- 
tion rate. Clearly it will be proportional to the surface 
area of the black hole, and to the energy density of the 
radiation background; the key question is the constant of 
proportionality, to which our results will prove extremely 
sensitive. The relevant length scales in the problem are 
the cosmological horizon size, the black hole radius, and 
potentially also the mean free path of the particles com- 
prising the radiation background (which should be much 
less than the horizon size if the background is to be ther- 
malized) . 

The two approaches available to estimate that constant 
are to use an absorption cross-section for radiation inci- 
dent from a distance, or to consider the thermal proper- 
ties of radiation near the event horizon. The former has 
the drawback that it assumes a long mean free path for 
the radiation (which cannot be true if the PBH radius is 
a significant fraction of the cosmological horizon), while 
in the latter the effect of the black hole geometry cannot 
readily be included. We will consider both and compare. 

To use the cross-section of the black hole, one needs to 
take into account that the impact parameter for absorp- 
tion of incoming radiation by the black hole is greater 
than the event horizon size, as the black hole bends the 
particle trajectory towards it. The black hole therefore 
has an effective radius r e g for capturing particles. In the 
standard cosmology r e g = 3\/3 r /2, where r = 2M/M% 
is the Schwarzschild radius ||. In the case of an effec- 
tively five-dimensional black hole, the equivalent expres- 
sion is r cS = 2r [Eol where 
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The accretion rate is estimated by assuming first of all 
that the radiation is non-interacting, so that the radia- 
tion reaching the black hole at time t has come from a 
known distance x. The fraction of the radiation starting 
at that point (assumed to have an isotropic momentum 
distribution) which is absorbed by the black hole is just 
the solid angle subtended by the black hole at that dis- 
tance, and then adding up over all particles in a shell of 
width dt at that distance yields 



dM 



(2) 



where the distance x drops out of the calculation. This 
expression will continue to hold even in the presence of 
interactions, at least as long as the particle mean free 
path is much larger than the black hole radius, because 
in thermal equilibrium particles are as likely to scatter 
onto trajectories towards the black hole as off those tra- 
jectories. However it is unclear precisely what effect scat- 
tering near the black hole might give, particularly in re- 
gions where the black hole has significantly modified the 
geometry, and so we can expect this result to be modified 
by a factor of order unity. It is also unclear whether the 



depletion of the radiation in the vicinity of the black hole 
by the absorption might lead to a significant reduction 
in accretion efficiency. Finally, we note that this calcu- 
lation is certainly naive in ignoring spin and frequency 
dependent effects in the absorption cross-section 0, 12| . 

A little more insight can be obtained by considering 
the thermal balance between the black hole and the back- 
ground if they were at the same temperature. The black 
hole radiation leads to a mass loss given by the Stefan- 
Boltzmann law applied to the radius r c g fic|| , except 
that it is also known that in the four-dimensional case 
the Stefan-Boltzmann law overestimates the emission by 
a factor 2.6 |TT[ , due to so-called grey-body factors al- 
lowing for the finite size of the black hole as seen by 
long-wavelength radiation (the majority of the power is 
emitted at wavelengths comparable to the event hori- 
zon radius) and the inclusion of spin-dependent effects. 
Such grey-body factors would also apply to absorption of 
long-wavelength radiation, and arguing that there should 
be balance when the cosmological and black hole tem- 
perature match suggests that in that situation Eq. (||) 
overestimates the absorption by a factor of around 2.6. 
However the regime of interest for accretion is when the 
background temperature is much larger than the black 
hole temperature, and then the grey-body factors should 
be much less important, supporting the normalization of 
Eq. (^). Were the cosmological temperature smaller the 
grey-body factors would be of increasing significance, but 
in that limit evaporation dominates in any case. 

Care is however needed in applying thermal balance ar- 
guments in the braneworld case, because while accretion 
is taking place on the brane, the black hole can evaporate 
into both brane and bulk; in effect the black hole provides 
a route for leakage of energy from the brane via evapo- 
ration of gravitons. Even if the black hole is in thermal 
balance on the brane, it will not be in equilibrium with 
the bulk, and so would still lose mass. However the ther- 
mal balance argument can be used between the accretion 
and the evaporation onto the brane alone. 

An alternative view is to consider the radiation back- 
ground at a given time to extend smoothly all the way 
to the event horizon with an isotropic momentum dis- 
tribution, and compute the flux entering the black hole. 
Consider a small shell of radius dt at the event horizon, 
which will contain an energy density AirrQp(t)dt. How- 
ever not all this radiation will be inwardly directed; the 
flux entering the black hole is one quarter of that energy 
density yielding 



dM o , . 

~dT = • 



(3) 



This is a factor of a few smaller than the estimate above. 
However this estimate is arguably more dubious. The 
approximation of an isotropic radiation distribution near 
the event horizon is unlikely to be good, as flux coming 
from the direction of the black hole will be absent, though 
one could argue that that this missing flux is correctly 
allowing for depletion of the radiation already absorbed 
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by the black hole. More seriously, this estimate does not 
allow for the effect of the black hole geometry on the 
radiation, whereas the effective radius for absorption did 
in the previous calculation. 

Summarizing the above, we can take the accretion rate 
to be 



(4) 



where F is a numerical constant measuring the accretion 
efficiency. It is conceivable that it might be as large as 
unity (particularly if the radiation mean free path is be- 
lieved to be much larger than the event horizon radius), 
whereas other arguments suggest it may be somewhat 
smaller. However neither calculation is likely to be ac- 
curate if the black hole is a significant fraction of the 
horizon size, because then accretion is constrained by the 
amount of material actually accessible to the black hole. 
Unfortunately, it turns out that the results in the high- 
energy regime are exponentially sensitive to the value of 
F, and therefore we keep it in the calculations that follow. 
By contrast, the recent paper by Majumdar ]4j assumed 
F = 1 throughout. 



III. PRIMORDIAL BLACK HOLE ACCRETION 
A. Accretion in the standard cosmology 

As a warm up, we can apply the above formalism in 
the standard cosmology, where however it will turn out 
that the approximations made are unlikely to be valid. 

Suppose a PBH is formed at a time U in the early 
phases of a radiation-dominated Universe, with mass 



M = fM H (U) 



fMlu . 



(5) 



Here Mh(U) is the horizon mass and / < 1 denotes what 
fraction of the horizon mass the black hole comprises. 
Using Eq. (Q) with the four-dimensional Schwarzschild 
solution and the radiation energy density 



P = 



3M 4 2 
327rf 2 



gives 



dM 
~dT 
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Integrating from a time t\ when the black hole mass is 
Mi leads to 



M(t) = Mi 



t/ti 



At/U + (1-A) 



A 



The behaviour of M(t) crucially depends on the sign 
and magnitude of the factor A. If it is positive, the black 
hole mass asymptotes to 



Mr, 



Mi 
~A 



(9) 



If F and/or / are small (small efficiency or initial PBH 
size), the factor will be close to 1, and the asymptotic 
mass will not be much bigger than the initial mass. But 
if A is small (Ff — > 0.4), the black hole mass will grow 
nearly as fast as the horizon mass, until tjt\ ss 1/A or un- 
til the radiation-dominated regime comes to an end. The 
factor A can even be negative for Ff > 0.4, in which case 
Eq. (||) gives a diverging mass in a finite time, violating 
causality and thus certainly indicating that the approxi- 
mations used have broken down. 

However, in order to obtain a final mass that is much 
bigger than the initial mass, it is required that the ini- 
tial PBH radius must be close to 74 percent (= 0.4 1 / 3 ) 
of the horizon radius (for efficiency F = 1), or more (for 
F < 1). For such large PBHs use of the quasi-static 
approximation (where the black hole is represented as a 
series of Schwarzschild solutions as the mass increases), 
becomes questionable. But in the case of PBHs the initial 
size of the black hole must be of the order of the hori- 
zon, to overcome the fluid pressure. Carr and Hawking 
Q proved there was no black hole solution embedded in 
a Friedmann background in which the black hole could 
grow at the same rate as the cosmic horizon. They then 
concluded that the black hole could only grow less fast 
than the horizon. After a short period of time its size 
would have become much smaller than the horizon, at 
which point the quasi-static approximation could be used 
to confirm no further accretion. This result has become 
widely accepted, though see Refs. (|, [l3| . 



B. Accretion in the high-energy regime of the 
braneworld scenario 



We now turn to the high-energy regime of the brane- 
world scenario. If black holes are to behave as five- 
dimensional ones, we showed [|J that they must form 
during the high-energy regime, and that it was possible 
for such PBHs to survive even to the present. In order to 
interpret possible observational signatures, we therefore 
need to understand accretion in the high-energy regime. 
In this section we consider only accretion, and in the fol- 
lowing section we will consider the combined effects of 
accretion and evaporation. 

If a black hole forms at a time i; in the high-energy 
regime of the Randall-Sundrum Type II scenario, its size 
will necessarily be smaller than the AdS radius I, and 
hence described as a five-dimensional Schwarzschild black 
hole. We take its initial mass to be a fraction / of the 
horizon mass (now computed using formulae relevant to 
the high-energy regime 0) 



M i =fM H (t l ) = WfM 4 



(10) 



The growth rate is again given by Eq. (^), where we must 
now employ the energy density in the high-energy phase 
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P = 



3 Mj 
32tt t c t ' 



(11) 



with the transition time between high-energy and stan- 
dard regimes given by t c — 1/2, and use the five- 
dimensional expressions for the event horizon radius and 
capture cross-section. This gives 



dM _ 2F M 

dt 7T t 



(12) 



M{t) = Mi 



2F/it 



(13) 



In this case we see that the black hole mass always grows 
less fast than the horizon mass M#(t) oc t 2 (provided 
F is not significantly greater than unity), in contrast to 
the situation in the standard cosmology. Accordingly, 
the approximations made in the calculation become more 
reliable as time goes by. Our result agrees with that of 
Majumdar who however assumed F = 1 throughout. 

This calculation shows that, regardless of the initial 
mass, black holes can experience significant growth, with 
the power-law depending sensitively on the accretion ef- 
ficiency F. The uncertainty in its precise value lends 
considerable uncertainty to the accretion growth of the 
PBHs during the high-energy phase. 

Once the high-energy phase ends, the cosmological 
background evolution becomes the standard one, but the 
PBH remains in the five-dimensional regime. To deter- 
mine whether accretion still continues, we now need to 
apply Eq. (Q) with the five-dimensional PBH properties 
but the standard cosmological evolution, taking as initial 
condition the mass of the PBH at the transition time t c 
and bearing in mind that the black hole will be much 
smaller than the horizon size at this point. Using these 
equations gives 



M(t) = M(t c ) exp 



2F 

7T 



1 - 



t, 



t 



(14) 



indicating a further growth by only a factor of order unity 
(comparable to the expected uncertainties in the calcu- 
lation) . 

We therefore conclude that accretion swiftly becomes 
unimportant as the standard cosmology sets in for five- 
dimensional black holes as well as four-dimensional ones. 
Accordingly, it is the slower decrease of the background 
density during the high-energy regime which makes ac- 
cretion important, rather than the change in black hole 
properties. 



IV. COMBINING ACCRETION AND 
EVAPORATION 

We now take both accretion and evaporation into ac- 
count in the high-energy regime 2 



where 



dM _ /dM\ ( dM 



(dM\ _ <? M 



cvap 



Integrating from the time of formation onwards results and 



(dM\ 



9 Ml 
2M 



(15) 



(16) 



(17) 



(see Ref. [||), where for convenience we have defined 

AF 



(18) 



Assuming the four-dimensional cosmological constant 
vanishes, the relation between the five-dimensional fun- 
damental mass scale and the Planck mass is given by 



Mr, = M i [ l - 



-1/3 



In the above expression we have defined 
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0.0062 

Cbranc 



0.0031 



9U 



Gbul! 



5bulk ■ 



(19) 



(20) 



where <?brane is the usual number of degrees of freedom 
into which the black hole can evaporate, while gbuik = 
0(1) is the number of bulk degrees of freedom (in the 
simplest case just the five polarization states of the gravi- 
ton). For the most part, the black hole's energy is lost 
through Hawking radiation on the brane. As an exam- 
ple we mention the case where the black hole emits only 
massless particles, for which gbranc = 7.25 and g = 0.023. 
Unlike in Ref. Q, we have written the grey-body factors 
Gbrane and Gbuik explicitly; in the standard cosmology 
the grey-body factor is equal to 2.6, but precise values are 
not known for the braneworld. For thermal balance be- 
tween the accretion and the evaporation onto the brane, 
the accretion efficiency F should equal 1/Gb ra nc when 
the temperatures are equal; however as discussed in Sec- 
tion II the absorption grey-body factors should approach 
one if the background temperature is much greater than 
the black hole temperature. 



2 Some work loosely related to this, but concerning astrophysical 
accretion by collisionally-produced mini black holes in a TeV 



gravity model, can be found in Ref. [Q. 
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In the standard cosmology Eq. ( [1 5[ ) does not have an 
analytical solution, but it does in the present case: 



M(t)/Mi = 




4v?l-9M ; 3/2 



(21) 

1/2 



An expression of this form was found by Majumdar 
though again for the specific case of efficiency F = 1 
(i.e. q = 4/71"). However, a very slight change in the ef- 
ficiency can change the qualitative behaviour. If q > 1 
(efficiency F > 0.78) the PBH will steadily grow until 
the cosmological transition time t c is reached (provided 
Mi > M 5 ). If q < 1 (efficiency F < 0.78), the mass loss 
term will grow faster than the mass gain term. However, 
assuming that the initial black hole mass is large com- 
pared to the fundamental mass scale, Mi ^> M5, the loss 
term will initially be much smaller than the gain term. 
In other words, if a black hole forms with a mass of order 
the horizon mass at a time t\ 3> tg , then the initial black 
hole temperature will be much lower than the temper- 
ature of the radiation background and evaporation can 
initially be neglected. 

For many choices of parameters, evaporation is then 
negligible throughout the high-energy regime. But for 
sufficiently low accretion efficiency, the mass growth can 
come to a halt whilst the black hole is still in the 
high-energy regime. 3 This 'halt' time is obtained from 
dM/dt = and reads 



! + (!-?) 



W7( 



Mi 



9 \M 5 



3/2' 



(22) 



If Mi S> M5 we can neglect the first term, to obtain 



1-9 



9(1-9) 



V7 /M 

9 \M 5 



3/2 



(23) 



In order for the derivation to be valid the halt time must 
satisfy t^ < t c , or equivalently 



4-g 



< 



9_ 

P 9(1 



2 g-9 



9) 



2-9 



(24) 



The above condition will be satisfied if I is large enough 
(a long high-energy regime) and/or if Mi is small enough 
(a high initial PBH temperature), depending on the value 
of q. Put another way, for a given I and Mi there is a 
minimum accretion efficiency which ensures that neglect- 
ing the mass loss through evaporation is justified all the 
way up to t = t c . 



3 In the limit q — » 1 the mass growth decreases logarithmically, 
and if Mj ^> M5 , evaporation can be neglected until the standard 
regime is reached. 



As an example we bear in mind that the AdS radius I 
is constrained by experiment as 



I < 10 31 L 



(25) 



If there was a period of high-energy inflation, the require- 
ment that gravitational waves do not lead to excessive 
anisotropies in the CMB leads to a lower limit on the 
horizon mass (and hence on the PBH mass) at the end 
of inflation, namely [p| 



Mj > 2 x 10 b M 5 



(26) 



For the extreme values of M; and I (and hence for all 
other values), we find that an efficiency better than 30% 
(q > 0.39) is sufficient to ensure mass growth all the way 
up to the cosmic transition time. 

We now ask what is the total lifetime t e vap of a black 
hole that reaches a halt time in the high-energy regime. 
Provided we can trust Eq. (^lj) until total evaporation 
M(t evap ) = 0, we find 



to 



(27) 



Unless q is very close to zero, we see that the total lifetime 
will be of the same order as the halt time (taking into 
account the possibility that the final stage of the black 
hole's lifetime is in the standard regime only decreases 
the estimate for i e vap> as it turns off the accretion term). 
Such black holes evaporate long before any observational 
constraints can be brought to bear, and we therefore con- 
clude that PBHs in the regime where accretion in the 
high-energy regime can halt are not of interest. 



To summarize, for a PBH with a lifetime t c 



^> t c we 



need consider only two situations in analyzing their evo- 
lution forwards in time. For black holes forming after the 
end of the high-energy regime, accretion is never impor- 
tant. For those forming within the high-energy regime, 
it is always a good approximation to neglect evaporation 
up until the cosmic transition time, after which accretion 
can be ignored and evaporation dominates. 



V. CONCLUSIONS 

We have explored the possibility that PBHs might ac- 
crete from the cosmological background in the case of 
braneworld cosmology, and found, in agreement with a 
recent paper by Majumdar [ij, that significant growth is 
possible. However, we have highlighted the extreme sen- 
sitivity of the resulting growth to the assumed accretion 
efficiency, which cannot be accurately computed. Accre- 
tion therefore adds considerable uncertainty to the evolu- 
tion of individual PBHs after formation in the braneworld 
scenario. 

Since accretion ends once the standard cosmology 
is restored (whether the PBHs are effectively four- 
dimensional or five-dimensional) , we stress that accretion 
does not have any implications for interpreting observa- 
tions in terms of the PBH density at evaporation. It does 
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however impact on how those constraints are interpret- 
ted in terms of formation rates in the early Universe. 
We will be providing a detailed analysis of observational 
constraints on braneworld PBHs in a forthcoming paper 
(Clancy, Guedens and Liddlc). 



Barrow and Malcolm Fairbairn for discussions. 
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